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Motivation 


•  Propulsion  systems  for  Nano  and  Pico-satellites 

-  0.1  to  10  kg 


-  High  energy  density  of  chemical  propellants  is  attractive 

-  Combustion  well  established  means  for  energy  release 


•  Behavior  of  deflagration  wave  may  be  different  when 
stabilized  in  small  passages 

-  Increased  exchange  with  structure 

•  Quenching 

•  Increased  burning  rate 

•  Enhanced  flame  stability 

-  Potentially  important  effects  on  performance  of  micro¬ 
rocket  motors 
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Previous  Work:  Fundamental 

UNIVERSITY  OF  MARYLAND 

•  Quenching  in  small  passages 

-  Zeldovich  (1941) 

-  Lewis  and  Von  Elbe  (1961) 

•  Effect  of  heat  loss  on  RR  and  flammability  limits 

-  Spalding  (1954) 

•  Heat  loss  decreases  burning  rate  and  broadens  flammability  limits 

•  Effect  of  heat  recirculation  (excess  enthalpy  burners) 

-  Weinberg  (1970) 

-  Weinberg  and  Hardesty  (1974) 

-  Takeno(l  979,81) 

•  Constant  T  walls,  single-step  reaction 

-  General  observations: 

•  Super-adiabatic  flame  temperatures 

•  Elevated  burning  rates 
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Previous  Work:  Micro-Channels 


UNIVERSITY  OF  MARYLAND 

•  Burning  rate  and  flammability  measurements  in  conductive  tubes 

-  Zamaschikov  (1997) 

•  Combustion  below  ‘quenching  limit’ possible 

•  Model  of  flame  propagation  in  a  narrow,  conductive,  channel 

-  Zamaschikov  and  Minaev  (2001) 

•  Fast  chemistry 

•  Conduction  broadens  flammability  limits 

•  Hysteresis  possible 

•  Minimum  ‘practical’  volume  for  an  HCCI  combustor 

-  Aichlmayr  (2002) 

•  Thermal  coupling  between  gas  and  structure 

•  Full  chemistry 

•  No  conduction  within  structure 

•  Effect  of  velocity,  heat  loss,  and  passage  width  on  burning  rate  in  microchannel 

-  Matalon  and  Daou  (2002) 

•  Constant  T  walls,  single-step  overall  reaction 

•  Heat  loss  to  environment  decreases  burning  rate 

•  Effect  of  axial  conduction  in  heat  recirculating  burner 

-  Ronney  (2003) 

•  PSR  with  full  chemistry 

•  Effect  of  axial  conduction  in  silicon  micro-combustor 

-  Current  work: 
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Obj  ecti  ve 


•  Investigate  physics  of  fluid  structure  coupling  that 
occurs  in  chemically  reacting  systems  operating  at 
micro-scales. 

-  Enable  development  of  small,  efficient  combustors  for 
micro  rocket  motors 

-  How  small  can  practical  micro-rockets  be  built? 

•  Efficient  and  stable  combustion 

-  What  level  of  performance  is  available? 

•  Thrust/weight  (power  density) 

•  Specific  impulse  (efficiency) 


6/21/04 


International  Symposium  on  Energy  Conversion  Fundamentals,  Istanbul 


UNIVERSITY 


Approach 


•  Modeling  and  Simulation 

-  Simple  model  for  effect  of  thermal  coupling  on  Sr 

•  Modification  of  Mallard-Le  Chatelier  approach 

-  Numerical  simulation 

•  1-D  geometry,  Full  chemistry,  Conjugate  heat  transfer,  Heat 
conduction  in  structure 

•  Experiments 

-  Develop  novel  non-intrusive  diagnostic  technique 

-  Investigate  behavior  in  parallel  plate  flow  reactor  with 
conductive,  temperature-controlled  walls 
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Analytical  model 


•  Include  thermal  coupling  with  structure: 


Structure  provides 
another  path  for 
heat  transfer 
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Analytical  model 


•  Thermo-electrical  analogy: 


Rs 
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Analytical  model 


•  Flame  thickness: 


with  P 


1+C 


k. 
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r  V 


S.  Nu 
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•  Asymptotic  behavior: 


-  Nu  — ►  0  or  H  — ►  00 :  No  thermal  coupling 


-  Nu  — > 00  or  H  — ►  0:  Perfect  thermal  coupling 


5r,fr  l  K 


Limit-cycle  behavior  of  5r  in  H 
Large  H:  limited  by  kj. 

(thin  flame) 

Small  H:  limited  by  ks 
(thick  flame) 
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Analytical  model 


•  Solving  for  8r  O  third  order  polynomial 


Questions: 
Influence  of  Nu? 
Influence  of  ks? 
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Numerical  Simulation 


T0=300K 
P=  1  atm 
H,-air 
<D=0.5 


Structure:  Energy 
Gas-phase: 

•  Continuity 

•  Energy 

•  Species 

Structure:  Energy 


201  grid  points 
Unevenly  spaced  mesh 
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VERS 


Numerical  Simulation 


•  Differences  wrt.  analytical  model: 

-  Chemistry  (9  species  and  1 9  reactions) 

-  Include  species  diffusion 

-  Transient  capability,  but  consider  steady-state  solution 

•  Assumptions: 

-  1-D 

-  P=const 

-  Laminar,  incompressible,  inviscid 

-  Nu=const 

-  No  thermal  diffusion 

-  Fourier  conduction  and  Fickian  diffusion 

-  Adiabatic  and  non-adiabatic  operation 
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Numerical  Simulation 


UNIVERSITY  OF  MARYLAND 
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Flame  Thickness 


•  Estimate  based  on  temperature  change  through  the  flame* 


*  Law,  C.K.  and  Sung,  C  .J.,  ‘Structure,  Aerodynamics,  and  Geometry  of  Premixed  Flamelets’, 
Progress  in  Energy  and  Combustion  Science  26  (2000)  459-505 
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Effect  of  Nusselt  Number 


•  Increasing  Nu  causes  broadening  to  occur  at  larger  H 

•  Consistent  with  model  predictions 
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Effect  of  Thermal  Conductivity 


•  Increasing  ks  increases  broadening  effect 

•  Consistent  with  model  predictions 
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Effect  of  Heat  Loss 


•  Heat  loss  reduces  broadening  and  leads  to  quenching 

•  Consistent  with  model  predictions 
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Silicon  Micro-Combustor 


VERS 


MARYLAND 


Conduction  also  increases  burning  rate. 
Including  k  JT)  is  important 
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Silicon  Micro-Combustor 


Definition  of  Power  Density 

Tt 

PSL\Cp(T)dT 


WD  = 


Note:  In  following  plots,  power  density  is  non-dimensionalized  by  a 
‘reference’  value  corresponding  to  combustion  at  the  laminar  flame 
speed  SL  in  a  volume  of  1  cm3. 
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Silicon  Micro-Combustor 


Power  density  increases  as  reduce  H;  optimum  L 
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Silicon  Micro-Combustor 


Definition  of  Efficiency 
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2/  H 
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Silicon  Micro-Combustor 


Adiabatic  Power  Density  and  Efficiency 


-1  0  1 

10  10  10 

L  /  8r,f  r  p=  1  atm 


Config.  That  maximizes  power  density  does  not  maximize  efficiency. 
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Continuum 


H/5 


Silicon  Micro-Combustor 

ERSITV  OF  MARYLAND 


Non-adiabatic  Power  Density  and  Efficiency 


10'1  10°  101 

L/5rfr  P=  1  atm 


Including  heat  losses  leads  to  optimum  H  and  L 
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Silicon  Micro-Combustor 


Non-adiabatic  Power  Density  and  Efficiency 


Increasing  pressure  increases  power  density 
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Silicon  Micro-Rocket 


VERSITY  OF  MARYLAND 
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Silicon  Micro-Rocket  Ln 

bp 


Maximum  I  at  large  H  and  large  L 
-  Reflects  trend  in  combustor  efficiency 


H/8 


yF  Silicon  Micro-Rocket  Thrust/Weight 


Practical  Design  Region? 


Optimum  Configuration 
(Neglecting  Pressure  Loss) 

•  L=0.5  mm 

•  H  =  5  pm 

•  T=72xl0-9  N 

•  T/W  =  943 


•  Maximum  T/W  corresponds  to  maximum  power  density 

•  Very  high  T/W  may  be  possible 
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Experiments 

UNIVERSITY  OF  MARYLAND  HUM 

•  Provide  data  to  validate  model  and  simulations 

Effect  of  Twall,  ksgruc!ure  on 

•  Burning  rate 

•  Reaction  zone  thickness 

•  Approach 

-  Construct  parallel  plate  reactor 

•  Conductive ,  temperature-controlled  walls 

•  Re-configurable  (0.5  mm  <  H  <  10  mm) 

•  Build  using  conventional  mfg.  processes  (avoid  MEMS) 

-  Develop  appropriate  diagnostic  techniques 

•  Measure  temperature  and  species  concentration 

•  Sub -millimeter  spatial  resolution 

•  Non-intrusive 
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Parallel  Plate  Reactor 
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MARYLAND 


Thermocouples 


Features: 

• Controllable  BC ’s 
• Detailed  Measurements 
•Temperature 
•Heat  Flux 


FTIR  interrogation  Volume 
(with  mask) 


'Species  cone. 
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Compute 
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Fuel  (CFLQHx)  Temperature  Controlled  Supports 


Traverse 
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►  x 


Efficiency, 

Temperature 


Efficiency 


►  x 
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Parallel  Plate  Reactor 
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Contraction 


Removable  section  : 
permits  study  of  developing 
or  fully  developed  flow 


Passage  heightt  ~7 
500  microns  -  5mm 


Cooling  passages 
Cartridge  heaters 

Plates 
Substrate 

SiC  layer 


Passage  length:  12  cm 


Passage  width:  7  cm 
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Parallel  Plat< 


Fuel  and  C02 

concentrations  n  Oxidizer 
Flame  thickness  4^ 


Fuel 
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FTIR  Proof  of  Concept 
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OF  MARYLAND 


Fourier  Transform  Infrared  Spectrometer  (FTIR) 


a  Traverse 
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1 — \ 
1 _ 
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J—  I 

LI 

i  ^ 

~2mm  — h  k- 


Si  wafers 

M 


Plenum 


Mask  increases  spatial  resolution 
FTIR  beam 


Interrogation  volume 


<=>  C3H8+Air 


Silicon  walled  micro -burner 
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FTIR  Proof-of-Concept 
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Silicon  Micro-burner  Operation 


Split  flame  attached 
to  Si  plates 
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Temperature  Calculation  from  CO 


1 .  Identify  peaks 


CO  spectra 
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Temperature  Calculation  from  CO 


UNIVERSITY  OF  MARYLAND 


2.  Fit  peaks 


Theoretical  fit  to  Experimental  Peaks  -  12m  Downstream,  phi  =  1.17 


Wave  Number  cm-1 


Er(vj) 


Nr={2j  +  l)- 


kT 


X(2/+i)? 

7=0 


Er(v,j) 


Er  (v,  j)  =  hc[Bj(j  +  \)-DJ2(j  + 1)2  ] 


kT 
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Temperature  (K) 


Axial  Temperature  Distribution 


Temperature  vs.  Downstream  distance 


Downstream  distance  (mm) 
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Present  Work 


Improved  Optics 

-  Better  FTIR 

-  Higher  throughput 
external  beam  path 

IR  Camera  for  Plate  T 

Improved  spectral 
interpretation 
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Present  Work 


Vertical  Traverse 


•  Improved  burner 

-  Adjustable  H 

-  Improved  traverse 

-  Improved  flow  control 
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Conclusions 

UNIVERSITY  OF  MARYLAND 

•  Structural  heat  conduction  has  important  effects  on 
performance  of  micro-combustors 

-  Increases  reaction  zone  thickness 

-  Increases  burning  rate 

-  Leads  to  optimum  power  density  configurations 

•  Rocket  motors  with  T/W  >  400  appear  possible 

•  Chemical  quenching  still  important 

•  Viability  of  micro-rockets  hinges  on  tradeoff 
between  T/W  and  Icn. 

•  Experimental  verification  ongoing 
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Future  Work 


•  Simulation 

-  Incorporate  radiation  boundary  condition 

-  Incorporate  surface  chemistry 

-  Investigate  different  fuels 

•  Experiments 

-  Measure  5r(H)  in  micro-burner 

-  Construct  parallel  plate  flow  reactor 

-  Replace  external  optics  with  optical  fiber  system 
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